

! 


JPL Publication 87-19 


err 






//J -w. <?<&. 

977f6~ 


X? 



Noiseless Coding for the 
Magnetometer 

Robert F. Rice Jun-Ji Lee 


(NASA-CR-18132 5) NOISELESS CODING FOB THE 
HAGNETOMETER (Jet Propulsion Lab.) 36 p 
Avail: N IIS HC AU3/HF A01 CSCL 22B 


G3/19 


N87- 29592 


Unclas 

0097495 




June 15, 1987 


NASA 

National Aeronautics and 
Space Administration 

Jet Propulsion Laboratory 

California Institute of Technology 
Pasadena, California 


JPL Publication 87-19 


Noiseless Coding for the 
Magnetometer 

Robert F. Rice Jun-Ji Lee 


June 15, 1987 


NASA 

National Aeronautics and 
Space Administration 

Jet Propulsion Laboratory 

California Institute of Technology 
Pasadena, California 


The research described in this publication was carried out by the Jet Propulsion 
Laboratory, California Institute of Technology, under a contract with the National 
Aeronautics and Space Administration. 

Reference herein to any specific commercial product, process, or service by trade 
name, trademark, manufacturer, or otherwise, does not constitute or imply its 
endorsement by the United States Government or the Jet Propulsion Laboratory, 
California Institute of Technology. 


ABSTRACT 


Future unmanned space missions will continue to seek a full understanding of 
magnetic fields throughout the solar system. Severely constrained data rates during 
certain portions of these missions could limit the possible science return. This publication 
investigates the application of universal noiseless coding techniques to more efficiently 
represent magnetometer data without any loss in data integrity. Performance results 
indicated that compression factors of 2:1 to 6:1 can be expected. Feasibility for general 
deep space applications was demonstrated by implementing a microprocessor breadboard 
coder/decoder using the Intel 8086 processor. 

The Comet Rendezvous Asteroid Flyby mission will incorporate these techniques in 
a buffer feedback, rate-controlled configuration. The characteristics of this system are 
discussed. 
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I. INTRODUCTION 


A full understanding of the magnetic fields throughout the solar system has been 
an important scientific objective since the inception of the deep space programnH3] t 
Both the Mars Observer (MO) and the Comet Rendezvous Asteroid Flyby (CRAF) missions 
are expected to fly sophisticated magnetometers to continue this investigation. 

The purpose of this publication is to provide a preliminary understanding of the 
potential advantages that may be derived from applying certain well-defined data com- 
pression techniques to magnetometer instruments. The work described herein was 
primarily performed in 1 985. The results were incorporated into the instrument proposal 
which was eventually selected for the CRAF mission^]. 

This publication is the second such study for science instruments. The results of 
a similar study for the Gamma Ray Spectrometer are described in Ref. 5. The latter results 
are currently being implemented for the MO mission. 

EXPERIMENT DATA GATHERING MODEL 

Figure 1 provides a simple interpretation of the magnetometer data gathering process. 
The magnetometer samples the three-dimensional magnetic field every r seconds along 
the spacecraft path of motion. 

Based on the characteristics of the ISEE-C Vector Helium Magnetometer described 
in Ref. 2, the measurement vector resulting at time tj can be denoted by 

Vj = r(i) x(i) y(i) z(i) (1) 

where the x(i), y(i), z(i) are m-bit representations of the three orthogonal components 
resulting from a linear quantizer operating over a range of signal values specified by 
the identifier r(i). 

The instrument described in Ref. 2 has the ability to manually or automatically 
switch between eight dynamic ranges, whereas the actual instrument built for MO or 
CRAF may have fewer ranges. In order to include the existence of r(i) within the overall 
coding problem, we let b r be the required length of r(i) in bits where 

1 < 0 r < 3 (2) 

depending on the number of instrument dynamic ranges. 


1 



z 



Fig. 1. Magnetometer Measurements. 


Then the baseline representation of the sampled magnetic field requires 

3m + P r bits/measurement vector (b/v) (3) 

The data rate resulting from continuous sampling is then 

3m + /8 

R, (m, /3 r , t) = [ b/s (4) 

r 


Desires and Constraints 

When instrument parameters have been sized to provide a basic maximum data- 
gathering capability compatible with the science objectives of both MO and CRAF , the 
following parameters result: 


m < 12 b/component 


(5) 
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I 


and 


1 /r < 20 vectors/s (v/s) 


( 6 ) 


Using Eq. (4) with 0 r = 3 we get a maximum instrument data rate of 


pmax w 800 b/s 


(7) 


Unfortunately, data rate constraints placed on the magnetometer instrument will 
be considerably lower than this. As a consequence, the parameters m and r would need 
to be altered significantly to meet this constraint. These basic trade-offs are graphed in 
Fig. 2. Certainly, such adjustments would likely reduce the science value of the returned 
measurements. 


MEASUREMENTS/s 



Fig. 2. Baseline Instrument Data Rates, Rj Versus m and r. 
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The remainder of this publication investigates data compression as a means of 
obtaining such required reductions in data rate with higher data integrity (i.e., less 
modifications to m and r). 

DATA COMPRESSION STRUCTURE 

The techniques previously applied to image data!®]"!®] an< j t ^e Gamma Ray Spec- 
trometer^] are equally suitable here. Figure 3 shows the general data compression 
structure. 
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Fig. 3. Overall Data Compression Structure. 


As shown, data compression generally includes two steps. The first, called “irre- 
versible preprocessing," consists of modification to the input data which cannot be undone 
by a decompressor. Such irreversible preprocessing could include such things as modifi- 
cations to sample quantization or effective sampling rate. Irreversible preprocessing can 
provide reductions in rate, but its consequence is a reduction in the fidelity of reconstructed 
data. 

On the other hand, the second step in Fig. 3, called "noiseless coding," can achieve 
a reduction in rate without reducing the fidelity of data encoded. In Fig. 3, this means 
that the output of the noiseless coder, C, can be "decompressed" to reconstruct B exactly 
(but not necessarily A). 

The output of the overall compression structure in Fig. 3 will therefore reconstruct 
with a quality corresponding to B, but operate at a rate corresponding to C. 
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Noiseless Compressor 


If we constrain 


B = A (9) 

we obtain a magnetometer compressor which results in exact reconstruction, that is, 
a "noiseless" magnetometer compressor. We will initially restrict attention to this form 
as indicated by the dashed line in Fig. 3. Any changes in m and r can be considered as 
changes to the instrument output parameters resulting in input A. 

Once we have defined and analyzed such a noiseless coder, we will return to the 
subject of irreversible preprocessing. An efficient noiseless coder will produce data rates 
which vary with data activity. For applications such as CRAF, which have fixed instru- 
ment data rate allocations, this necessitates some form of "rate control." We will 
describe a buffer feedback system for CRAF which accomplishes this objective by mak- 
ing fine adjustments to the effective sampling rate (i.e., varying r)^. 


II. NOISELESS MAGNETOMETER COMPRESSOR 


The noiseless coding required here, as for the Gamma Ray Spectrometer^ 5 !, is a 
direct application of techniques described in Refs. 6 through 9. Whenever possible we 
will use the notations and definitions already described in these reports. 

PRELIMINARIES 

Notation 

Key notational quirks are worth noting before proceeding: 

1 ) We will generally place a tilde ( ~ ) over a symbol to indicate a sequence of bits 
or other symbols; 

2) An asterisk (*) will be used to indicate concatenation of symbols or sequences 
when we need some additional emphasis; 

3) If Z is a sequence, we use £( Z) to indicate its length in bits; and 

4) We subscript and superscript ^[»] to indicate a particular noiseless code 
operator. Specifically, if Z is an input sequence to some code operator ^J*], 
then \p [ Z] is the resulting coded sequence and, given \J/J Z], the original sequence 
Z can be recovered precisely. 


General Structure 

The noiseless coder in Fig. 3 generally splits into separate operations as indicated 
in Fig. 4. 

In addition to formatting incoming data, the function of reversible preprocessing 
is essentially to predict sample values and then map the prediction error into the integers 
such that the most likely occurring integers are the smaller ones^U^- The better this 
last objective is achieved, the better the variable length coder that follows will work. 
That is, a better predictor will lead to lower average coded rates. The process is revers- 
ible in that the original data can be recovered from the mapped integers. 
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Fig. 4. Overall Noiseless Coding Structure. 


The step of reversible preprocessing is generally source dependent, whereas the 
step of variable length coding can generally be made source independent. Its implemen- 
tation usually consists of selecting the appropriate algorithms from Refs. 6 and 7, 
which efficiently code integers over the range of data activity expected for the given 
application. 


SEPARATING THE DATA SOURCE 

Magnetometer data as defined in Eq. (1 ) really consists of four separate data sources 
which need to be dealt with individually. It will be convenient later to have isolated this 
step of preprocessing, and hence the following discussion will lead to an intermediate 
definition of a magnetometer compressor. 

Referring to Fig. 1, let 


V = v 1 * v 2 * . . . * Vj (10) 

be a block of J consecutive measurement vectors starting at time tj. That is, using Eq. (1 ) 

V = r( 1 )x( 1 )y ( 1 )z( 1 ) * r(2)x(2)y(2)z(2) 

* r(J)x(J)y(J)z(J) (11) 
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Now collect all the range and the x, y, z components into individual blocks 

R = r(1) r(2) . . . r(J) 

X = x(1 ) x(2) . . . x(J) 

( 12 ) 

Y = y(1) y(2) . . . y(J) 

Z = z(1) z(2) . . . z(J) 

and concatenate them into the new sequence 

V' = R * X * Y * Z (13) 

Clearly these rearrangements are reversible; V can be recovered from V'. The steps 
are illustrated in Fig. 5. 




z 


Fig. 5. Separating Into Multiple Data Sources. 


Note that V' has the sarr-e form as the original single measurement vector Vj in Eq. 
(1), except that now each segment is a J sample sequence of the same kind of data. 

We can now treat local segments of the data sources for range and x, y, z by 
separately coding the blocks R, X, Y and Z at point ® in Fig. 5. Using the code oper- 
ator notation of Refs. 6 and 7, let 

* R W. * x H. lAy [ • ] r (14) 

represent the noiseless code operators to be applied to R, X, Y and Z, respectively. 

From Fig. 5 and Eq. (14), the magnetometer noiseless coder (denoted ^|y|[*]) takes 
the form shown in Fig. 6. We need only identify the operations of the various ^[*]s. 
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Fig. 6. Magnetometer Noiseless Code Operator, ^/jM- 


The output of the magnetometer code operator as shown in Fig. 6 is defined by 

* M [V] = 0 r [R] * iA x [X] * * y [Y] * i z [Z] (15) 

and, since each of the individual operations are reversible (by definition), R, X, Y and Z 
can be recovered from ^[V], from which V itself can be retrieved. 

CODING X, Y AND Z 

Although X, Y and Z need to be treated separately, they can each use the same 
coding process, which at this point is identical to an image compressor as defined in 
Ref. 7. For convenience, let 


(16) 

be the general purpose code operator which applies to X, Y and Z so that 

* X M - ^ y [*3 ^ H7) 

Now let A be an input sequence to (i.e., X, Y or Z), and let a(j) be the j** 1 
element of A so that 


A = a(1 ) a(2) . . . a(J) 


(18) 
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Now if we use a(j- 1) as the prediction of <*j we get the j** 1 prediction errort 

A. = a(j) - a(j — 1 ) (19) 

and a J sample block of prediction errors 

A = A 1 A 2 . . . Aj (20) 

If we then map each of the individual errors, Aj, into the integers we get 

8 = 5 1 b 2 . . . Sj ( 21 ) 

where 

Aj - 8 . (22) 

according to the relationship in Table 1. 


Table 1. Basic Mapping of A Into the Integers 5. 


PREDICTION 


ERROR 

INTEGER 

A 

8 

0 

0 

+ 1 

1 

-1 

2 

+ 2 

3 

-2 

4 

+ 3 

5 

-3 

• 

6 

• 

• 

• 

• 

• 


tThe last sample of the previous A block can be used for a(0). If a previous block does 
not exist, a{0) can be a constant. 
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The latter steps are precisely the same steps shown in more detail in Ref. 7 and 
constitute the "reversible preprocessing" of code operator The remaining step 

is provided by a variable length coder which we will denote here as 

[•] (23) 

and which can be chosen from numerous alternatives in Refs. 6 and 7. 

The code structure for handling X, Y or Z through ^\[«] is shown in Fig. 7. 



H 





Fig. 7. Code Operator for X, Y, or Z. 


Modifications for Frequent Range Changes 

Note that when a range change occurs, the samples a(j - 1 ) and a(j) are generated 
from different dynamic ranges. The error Aj in Eq. (19) could be quite large even if the 
corresponding analog sample values were the same. The net effect of a large Aj (a spike) 
would be an increase in code rate for the block in which it occurred. All our test data 
showed a negligible effect on long-known average coder performance. 

However, in the event that more frequent range changes became the norm, predic- 
tion errors could be significantly reduced by mapping the prediction value a(j — 1 ) from 
its range into the range of the next sample. Further, Ref. 7 provides supplemental 
techniques for improving efficiency in data blocks containing spikes (e.g., occasional 
large errors). 
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CODER FOR RANGE 


The same principles apply to range data, although its characteristics are quite 
different than that of the x, y and z components. Changes in range are quite infrequent, 
so that one would expect that the rate required to communicate these changes should 
be considerably smaller than the 1 to 3 b/sample assumed in the basic representation of 
Eq. (2). This is in fact the case. 


Our primary test data had eight ranges, so we will make that assumption here 
although actual CRAF or MO instruments will probably have fewer ranges. The code 
operator defined below for range data can be viewed as simply an expedient definition 
which, from our limited test set, achieves average rates which are negligible compared 
to that of X, Y and Z. We have also built in some basic protection against unexpected 
conditions. Additional sophistication to improve efficiency would certainly be unwarranted. 


Let 


and 


f(i) 


0 if r(i) = r(i — 1 ) 

1 otherwise 


(24) 


w(i) = 


</> if f(i) = 0 
r(i) otherwise 


(25) 


where </>= non-existent, null or “no data" (see below). Then define code operator by 


Iro 


V'pt R] - r(1) co(1) f(2) co(2) . . . f(J) co(J) 


(26) 


In words, the £*(j) are 0's and 1 's indicating if a range has changed, and if a change has 
occurred (f = 1 ) that f is followed by a three-bit (in general /3 r ) specification of the new 
range. If a particular f equals zero, a new range value is not needed (co = </>) and that 
f is instead followed by the next f. 

1 

Code operator i^pH has two areas of inadequate performance. In the first, if all 
the range samples changed, \J/ p[«] would require four b/sample. Although unlikely, it is 
worthwhile to protect against this transient situation. 


First recall that by the notation of Refs. 6 and 7, and earlier discussions, if K is a 
sequence 


£( K) 

tRef. 7 also describes a slightly more efficient mapping algorithm. 


( 27 ) 
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means the length of that sequence in bits. Now define 


and let 


e = 


1 if £ ^ [ R] > 3J 
0 otherwise 


*r[R] = 


0 * R] if e = 0 

1 * R if e = 1 


( 28 ) 


(29) 


o ~ i ~ 

A close look at this operator reveals that ^r[ R] only uses ^p[»] to code R when it 
actually provides a gain over no coding at all. The worst-case code rate for ^p[«] is 
(3 + 1/J) b/sample. 

1 2 

The second deficiency in (and the one remaining in ^p[«]) is that its minimum 
rate of one b/sample is too high. This occurs in the very likely event (by our test data) 
that no range changes are detected in a J sample block. 

Now let 


f = 


and define the operator ^p[ # ] by 


<$R] = 


0 if no range changes 
in complete block 

1 otherwise 


0 if / = 0 

1 * *r[R] if / = 1 


which by Eqs. (28) and (29) expands to 


(30) 


(31) 


*r[R] 


0 if / = 0 (no changes) 


1 * 


0 * ^p[ R] if e = 0 

1 * R if e = 1 


(32) 
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3 ~ 3 ~ 

Looking closely at Eq. (32) we see that the length of ^r[R], £(^r[R]) is only one 
bit per block if there are no range changes (a high probability); and 2 + J + 3 (number 
of changes) up to a maximum of (2 + 3J) bits when the number of changes is greater 
than zero. 

This is illustrated in Fig. 8 for a block size of J = 16. 

OVERALL NOISELESS CODER 

We can now expand the magnetometer noiseless code operator, ^[«] in Fig. 6, into 
the more complete description in Fig. 9. 
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III. DATA COMPRESSION TESTS 

ISEE-C DATA 

The first magnetometer test set available to us was actual decalibrated data from 
the ISEE-C experiment described in Ref. 2. The data were represented by 

i 

I 

m = 9 b/component (33) 

which could be applied to eight dynamic ranges so that 

0 r = 3 (34) 


The sampling rate was at 

1/r = 6 vectors/s (v/s) (35) 

which is roughly a factor of three slower than the desired 20 v/s, as noted in Eq. (6). 
While not a perfect match to the anticipated configurations for CRAF and MO, the test 
runs clearly provided a feel for first-order effects with some reasonable extrapolations 
possible. 

Variable Length Coders 

The only undefined element of the block diagram in Fig. 9 is the variable length 
coder, designated ^ [•]. The three algorithms tested below are described in Refs. 6 
and 7. These are: 

A) An eight-option version of the "Fast Compressor" used at the recent 
Uranus encounter for Voyager imaging. 

B) A four-option adaptive coder called the "Basic Compressor" or [•]. 

C) An eight-option ^-| -j [•] using (see, for example. Fig. 10 in Ref. 

7). This version, while still simple and efficient, provides enough dynamic 
range to satisfy the requirements of most magnetometer instruments 
as well. 

Since these algorithms are essentially source independent and thoroughly discussed 
in Refs. 6 and 7, we will not provide any more details here. 
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Performance Measurements 


Average measured performance of some noiseless code operator will mean the 
average number of bits per (x, y or z) component required by that operator over some 
specified measurement span (i.e., sum the bits and divide by the number of measure- 
ments). We will look at such average performance from several different angles. 

Average entropy calculated for some measurement span is a practical guide to the 
minimum possible rate one can expect to code the corresponding data and still remain 
noiseless. Under some conditions it is an absolute bound to performance [6], [7] |-| eref 
entropy calculations (-E Pj log pj) were derived from distributions of prediction errors, 
A in Eq. (20) and range samples r(i). 

Overall Average. We first look at the average code rates obtained over the complete 
test set by operator ^[«] using the variable length code options A, B and C and ^p[«] 
as specified in Eqs. (24) through (32). The results are given in Table 2. 


Table 2. Overall Average Performance of ^[»]. 


ENTROPY 

(BITS/COMPONENT) 

MEASURED PERFORMANCE 
(BITS/COMPONENT) 

= OPTION A 

^[•] = OPTION B 

i P e ['] = OPTION C 

3.56 

3.88 

3.45 

3.42 

INCLUDES A 
RANGE ENTROPY 
OF « 0.06 BITS 
PER RANGE 
SAMPLE 

INCLUDES AN AVERAGE CODE RATE OF 0.11 BITS PER 
RANGE SAMPLE USING in EQS. (24) TO (32). 
















i We see that Option C reduces the fixed rate of 1 0 b/component (30 b/v) to only 

3.42 b/component on the average. This is an average compression factor of 2.92. 

i 

Variability. As noted earlier, the code rate actually achieved with an efficient noise- 
less coder will vary as the activity of the data (entropy) varies. The extent of this varia- 
tion is important to buffering requirements and/or feedback control of quantization and 
sampling rate to meet data system rate limitations. We take a first order look at this 
phenomenon by comparing the entropies for X, Y and Z components and their com- 
posite averages taken over 10 test blocks of 44,000 vector samples each. The result 
is plotted in Fig. 10. 

As the diagram shows, X, Y and Z activities follow the same trend in variation. 
In each case the observed swing between maximum and minimum entropies is about 
2 b/component. The same fluctuations are thus observed on the composite average of 
the X, Y and Z graphs. 

i 



x 

z — 

COMPOSITE - — 


Fig. 10. Rate Variability Over Ten 44,000 Vector Test Blocks. 
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VENUS IONOSPHERE DATAnO] 


Chris Russell of the University of California at Los Angeles provided us with several 
files of magnetometer data for the Venus ionosphere. The files contained data quantized 
to 12 b/component at sampling rates up to 24 v/s. 

Applying the same noiseless coding techniques, we observed average compression 
factors of 3:1 . Compression factors as high as 6:1 were observed over short stretches. 

GIACOBINI-ZINNER COMET DATA^l 

Further investigations included the use of 1 985 observations taken in the turbulent 
tail regions of Comet Giacobini-Zinner. These tests were used to evaluate a buffer- 
feedback, rate control system for the CRAF mission. The results are discussed in the 
next chapter. 
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IV. CRAF - A PRACTICAL APPLICATION 


The noiseless coding techniques just discussed formed a part of the magnetometer 
instrument proposal which was accepted for inclusion in the CRAF science payloadt^O. 
The following discussion looks at this practical application in more detail. 

AN OPEN LOOP COMPRESSION SYSTEM 

An expanded magnetometer compression system block diagram and some additional 
notation are shown in Fig. 1 1 . 

Raw Instrument Data Rates 

The proposed instrument will generate 20 vectors/s (v/s) using 1 2-bit quantization 
on X, Y and Z components. This and additional range information results in a raw instru- 
ment data rate oft 


Rj * 800 b/s (36) 

as in Eq. (7) for the baseline system. 

Filtered Rates 

A filter will be included which can provide selective adjustments in vector sampling 
rates by the factors 


/ s = 1 , 2, 3, 4 and 6 (37) 

yielding effective vector sampling rates of 20, 10, 6.67, 5 and 3.3 v/s, respectively, and 
corresponding filter output rates of 

R F = 800, 400, 267, 200 and 1 33 b/s (38) 

The inclusion of this filtering can be interpreted as a form of "irreversible preprocessing" 
in Fig. 3 or simply modifications to the instrument output characteristics. 


t Adjusted slightly for convenience in interpretation. 
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Code Operator, 

The noiseless coder from Fig. 9 (using ^[«] m Option C in Table 2) will be applied 
to the output of the sampling filter. We denote 

R c b/s (39) 

as the output rate of the compressor and 


e 


as the compression factor 



(40) 


(41) 


But 0 is also 


Uncompressed 

0 = ^/v (42) 

Compressed 

b/v 

which we know will vary over time as data characteristics vary. But for the same seg- 
ment of data, a change in sampling rate will not significantly alter the per vector perfor- 
mance of Thus 0 will remain roughly the same for each sampling rate (when applied 

to the same segment of data). Rewriting Eq. (41) we see that 


R 



(43) 


will vary proportionately with changes in filter output rate Rp. 

For example, if ^[*] generates Rq = 300 b/s when Rp = 800 b/s (20 v/s), it will 
generate Rq * 1 50 b/s if the filter rate is reduced by a factor of two to Rp = 400 b/s 
(10 v/s). 


Link Rates 


The communication link rate 


R L b/s 


(44) 
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which is made available to the magnetometer will not be fixed. In addition to the primary 
rate of 400 b/s, the magnetometer will be forced to operate with as little as 200 and 
100 b/s during various phases of the mission. 

Buffer 

A "first-in first-out" (FIFO) buffer will be included to handle: 


1 ) The variability in noiseless coder output rate, Rq; and 

2) The variability and limitations imposed by changes in link rates, Rj_. 


The FIFO buffer accumulates data from the noiseless coder. Simultaneously, it also 
continually transfers the oldest of this accumulated data to the data link at a rate of R|_ b/s. 

If the average compressed rate, Rq, exceeds the available link rate, R|_, the buffer 
will gradually fill up. If it did fill up, subsequent data from would be lost until enough 

room was provided by transfers of old data. The complete loss of data in this manner 
is viewed as quite unacceptable. 

But note that if Rq < Rl the buffer fullness will decrease. Thus, we should be able 
to avoid buffer overflow by adjusting the filter sampling rate downward to lower Rp and 
hence Rq = Rp/0 as needed. The latter approach is discussed next. 

CLOSED LOOP COMPRESSION SYSTEM 


The block diagram of Fig. 1 1 has been modified to include a feedback loop to the 
filter, controlled by a measure of "buffer fullness." The result is shown below in Fig. 1 2. 


RATE-CONTROLLED COMPRESSOR 
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Although a precise algorithm has not yet been specified, the basic idea is to switch 
the filter to a lower sampling rate (lowering Rp) whenever the buffer might overflow. 
Conversely, the sampling rate can be increased (higher Rp) when the buffer fullness gets 
too low. The fullness level itself, as well as the rate of change in buffer fullness, can 
be used as control parameters. 


Average Effect 

The long-term statistical effect of this control strategy is to maximize the average 
sampling rate. Lets look at a concrete example. 

Suppose the available link rate is R|_ = 200 b/s where an uncompressed data 
stream would require the sampling rate to be reduced to only 5 v/s. If the compression 
factor achieved were 4:1, a constant buffer level could be maintained at the full 20 v/s 
while matching the link rate R|_ = 200 b/s. 

If the compression factor were only 2: 1 , a constant buffer level could be maintained 
only if the sampling rate were reduced to 10 v/s. At any inbetween compression factor, 
the feedback mechanism would automatically cause switching between sampling rates 
such that the overall average compressed rate Rq (which is also the buffer input rate) 
equals the available link rate R|_, 200 b/s (which is the buffer output rate). 

At a 3:1 compression factor, about half of the transmitted data would be at 20 v/s 
and half would be at 10 v/s. 

Tests of this form were performed using 1985 observation data of the turbulent 
tail regions of Comet Giacobini-Zinnerf^]. Similar average results with stable buffer levels 
were observed. 

Transient Advantages 

As these tests indicated, it is rare that the compression factor achieved actually 
dips below 2:1 . At a link rate of 400 b/s, a constant buffer level can be maintained at 
the full 20 v/s sampling rate. 

A number of interesting magnetometer events can produce brief bursts of very active 
data with corresponding compression factors dipping below 2: 1 Q Ut jf the buffer 
maintenance level under normal conditions (0 > 2) is kept fairly low, the remainder of 
the buffer can be used to absorb a good part of such a burst at the full 20 v/s sampling rate. 
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More complex strategies would even allow portions of such bursts to be preserved 
at high sampling rates at the lower link data rates. 

OTHER MODES 

The discussion above describes the principal CRAF magnetometer application of 
noiseless coding. Other modes of operation will be discussed in later reports. 
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V. BREADBOARD 


A breadboard coder/decoder of in Fig. 9 (using ^[«], Option C in Table 2) 

was built and tested by Bruce Parham using an Intel 8086 microprocessor^ 1 The bread- 
board coder required roughly 1 000 bytes of memory for instructions and another 1 000 
bytes for internal buffers. A decoder was about half the size. Both handled throughput 
data rates of 20 kb/s, far higher than the allowable rates on CRAF or MO. 
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